## Drawing Fig9(a)
import matplotlib.pyplot as plt
from matplotlib import rcParams

plt.rcParams['figure.dpi'] = 1200
fr = open("DrawingFig9(a)(b).csv", "r")
ls = []
for line in fr:
    line = line.replace("\n","")
    ls.append(line.split(","))
fr.close()

x_start = 0.0
x_end = 0.20
y_start = 0.0
y_end = 0.07
x = []
y = []
iLen, jLen = len(ls), len(ls[0])
for i in range(1, iLen):
    x.append(float(ls[i][2]))
    y.append(float(ls[i][4]))

plt.figure(figsize=(6.0, 4.0))
plt.xlim(x_start-0.01, x_end+0.01)
plt.ylim(y_start-0.005, y_end+0.005)
plt.xticks([0.00, 0.05, 0.10, 0.15, 0.20], ['0.00', '0.05', '0.10', '0.15', '0.20'])
plt.plot(x, y, "g.-", linewidth=1)
plt.xlabel('$r/L$')
plt.ylabel('$Std$')

plt.savefig("Fig9(a).jpg")
plt.show()
